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■ We relate the graph isomorphism problem to the solvability of certain sys- 

\ terns of linear equations and linear inequalities. The number of these equations 

and inequalities is related to the complexity of the graphs isomorphism and 
subgraph isomorphim problems. 

2000 Mathematics Subject Classification: 03D15, 05C50, 05C60, 15A48, 15A51, 
15A69, 90C05. 

| Keywords and phrases: graph isomorphism, subgraph isomorphism, tensor 

^ ■ products, doubly stochastic matrices, ellipsoidal algorithm. 

1 Introduction 

m 
> 

Let G\ = (V, E\),G-i = (V, E%) be two simple undirected graphs, where V is the set 
of vertices of cardinality n and E\,Ei C VxV the set of edges. G\ and Gi are called 
isomorphic if there exists a bijection a : V — > V which induces the corresponding 
bijection a : E\ — > E^. The graph isomorphism problem, abbreviated here as GIF, 
is the problem of determination if G\ and G2 are isomorphic. Clearly the GIP in the 
class NP. It is one of a very small number of problems whose complexity is unknown 
[4, 6j. For certain graphs it is known that the complexity of GIP is polynomial 

m ei s ei m Eg. 

^ . Let G3 = (W,Es), where #W = m <n. G3 is called isomorphic to a subgraph 

of G2 if there exits an injection r : V3 — > V% which induces an injection f : E3 — > Ei- 
The subgraph isomorphism, abbreviated here as SGIP, is the problem of determi- 
nation if G3 is isomorphic to a subgraph of G%. It is well known that SGIP is 
NP- Complete [3]. 

In the previous versions of this paper we related the graph isomorphism problem 
to the solvability of certain systems of linear equations and linear inequalities. It was 
pointed out to me by N. Alon and L. Babai, that my approach relates in a similar 
way the SGIP to the solvability of certain systems of linear equations and linear 
inequalities. Hence f(n), the number of these linear equalities and inequalities 
for V = n, is probably exponential in n. Thus, the suggested approach in this 
paper does not seem to be the right approach to determine the complexity of the 
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GIP. Nevertheless, in this paper we summarize the main ideas and results of this 
approach. It seems that our approach is related to the ideas and results discussed 
in [IT]- 

Let n C IR™ X?1 be the convex set of n x n doubly stochastic matrices. In this 
paper we relate the complexity of the GIP to the minimal number of supporting 
hyperplanes determining a certain convex polytope ^> n ,n C £l n 2 . More precisely, two 
graph are isomorphic if certain system of n 2 hyperplanes intersect \P nn . More gen- 
eral, if the corresponding system n 2 half spaces intersect ty n n then G3 is isomorphic 
to a subgraph of G<i- Hence the minimal number of supporting hyperplanes defining 
^ n ,ni denoted by f(n), is closely related to the complexity of SGI P. We give a 
larger polytope $ n ,n> characterized by {An — l)n 2 linear equations in n nonnegative 
variables satisfying 

^n,n C $ n , n C Q n 2. (1.1) 

In the first version of this paper we erroneously claimed that & n ,n = ^ n ,n- The 
error in my proof was pointed out to me by Babai, Melkebeek, Rosenberg and 
Vavasis. The inequality ^f n ,n ^ ^n,n for ra > 4 is implied by the example of J. 
Rosenberg. 

Thus if two graphs are isomorphic then certain system of n 2 hyperplanes intersect 
$ n ,n- This of course yields a necessary conditions for GIP and SGIP. 

We now outline the main ideas of the paper. Let A,B be n x n adjacency 
matrices of G\,G2- So A, B are — 1 symmetric matrices with zero diagonal. It is 
enough to consider the case where A and B have the same number of l's. Let V n be 
the set ofnxn permutation matrices. Then G\ and G2 are isomorphic if and only 
if PAP T = B for some P G V n . It is easy to see that this condition is equivalent 

P(A + 2I n )Q T = B + 2I n for some P, Q G V n , (1.2) 

where I n is the nxn identity matrix. 

For C,D G R nxn denote by C ® D G R™ 2 *™ 2 the Kronecker product, see §2. 
Let V n <g> V n := {P <g> Q, P,Q G V n }- Denote by ^ n ,n C M+ x ™ 2 the convex 
set spanned by V n <S> V n . *& n ,n is a subset of n 2 x n 2 doubly stochastic matrices. 
Then the condition (jl.2|) implies the solvability of the system of n 2 equations of 
the form Z(A + 2I n ) = B + 2I n for some Z G ^ n , n . Here B + 2I n G M™ 2 is a 
column vector composed of the columns of B + 2I n . Vice versa, the solvability of 
Z{A + 2J n ) = B + 2I n for some Z G ^? n ,n implies fll .2[> . The ellipsoid algorithm 
in linear programming [HI [7j yields that the existence a solution to this system 
of equations is determined in polynomial time in max(/(n),n). Similarly, for the 

SGIP one needs to consider the the solvability of Z(C + 2 n2 I n ) < B + 2 n2 I n for some 
Z G ty n ,n, where C is the adjacency matrix of the graph G3 = (V, E3) obtained from 
G3 by appending n — m isolated vertices. 

We now survey briefly the contents of this paper. In §2 we introduce the needed 
concepts from linear algebra to give the characterization of <3? njn in terms of (4n— 2)n 2 
linear equations in n 4 nonnegative variables. This is done for the general set $ mn , 
which contains \P TOj n, the convex hull of V m ®V n - §3 discusses the permutational 
similarity of A, B G M nxn and permutational equivalence of A, B G M nxm . We 
show the second main result that the permutational similarity and equivalence is 
equivalent to solvability of the corresponding system of equations discussed above. 
In §4 we deduce the complexity results claimed in this paper. 
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This paper generated a lot of interest. I would like to thank all the people who 
sent their comments to me. 



2 Tensor products of doubly stochastic matrices 

For m G N denote (m) := {1, . . . , m}. For ^Cl denote by A mxn the set ofmxn 
matrices A = [dij]^-! such that each Oij G A. Recall that A = [a^] G M™ xm is 
called doubly stochastic if 

m m 

^Z a ij = ^Z a ji = 1j i = l,...,m. (2.1) 
i=i i=i 

Since the sum of all rows of ^4 is equal to the sum of all columns of A it follows that at 
most 2m — 1 of above equations are linearly independent. It is well known that any 
2m — 1 of the above equation are linearly independent. Let 1 := (1, . . . , 1) T G W+. 
Note that A = [ay] G M mxm satisfies fl23) if and only if and Al = A T 1 = 1. 

Denote by f2 m the set of doubly stochastic matrices. Clearly, Q m is a convex 
compact set. Birkhoff theorem claims that the set of the extreme points of O m is 
the set of permutations matrices V m C {0, l} mxm . 

Lemma 2.1 Denote by A m C M™ xm the set of nonnegative matrices satisfying 
the conditions Al = A T 1 = al for some a > depending on A. Then A m is a 
multiplicative cone: 

A m + A m = A m , aA m C A m for all a > 0, A m ■ A m = A m . 

Furthermore, A = [aij] G M™ xm is in A m if and only if the following 2(m — 1) 
equalities hold. 

m m m 

^ a-ij = ^2 a i i = ^2 ai i for i = 2,. . . ,m. (2.2) 
j=i j=i j=i 

Proof. The fact that A m is a cone is straightforward. Since I m G A m we 
deduce the equality A m ■ A m = A m . Observe next that the conditions (|2.2p imply 
that Al = al, where a is the sum of the elements in the first row. Also the sum 
of the elements in each column except the first is equal to a. Since the sum of all 
elements of A is ma it follows that the sum of the elements in the first column is 
also a, i.e = al. □ 

For A = [oij] G W mxm , B = [b M ] G R nxn denote by A B G M. mnxmn the tensor 
product of A and B. The rows and columns of A <g) B are indexed by double indices 
{i, k) and (j, /), where i,j = l,...m, k, I = 1, . . . , n. Thus 

A®B = [c( !l% i)]eR m " xmn . (2-3) 

where C(i,k)(j,i) = <Hjhl for i, j = 1, . . . ,m, k, I = 1, . . . , n. 

If we arrange the indices (i,k) in the lexicographic order then A ® B has the 
following block matrix form called the Kronecker product 
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A® B 



CL\\B CL12B 

a<i\B (122B 
a m \B a m 2B 



a\ m B 

&2mB 
Qm.m.B 



(2.4) 



For simplicity of the exposition we will identify A® B with the block matrix ([27 
unless stated otherwise. Note that any other ordering of (m) x (n) induces a dif- 
ferent representation of A (g) B as C G M mnxmn ) where C = P(A <g) B)P T for some 
permutation matrix P G V m n- 

Recall that A <g> B is bilinear in A and -B. Furthermore 



(A ® £)(C ® D) = (AC) (8 (BD) for all A,Ce M mxm , B, D G 
Proposition 2.2 Lei A G Q TO , B G 0„. T/ien i®£ 6 Q mn . 



(2.5) 



Proof. Clearly A ® is a nonnegative matrix. Assume the representation 
(1231). Then 



to n 



m,ra m,ri 

52 c (i,fc)0',0 = 0i i &fc * = (52 a v)(52 = 1-1 = 1, 

j,l=l j,l=l j=l 1=1 

m,n m,n m n 

52 c uM,k) = 52 a i ihik = (52 a ^(52 bik ^ = 1-1 = 1. 

3,1=1 3,1=1 j=i 1=1 



□ 



Lemma 2.3 Denote by ^ m ,n C £l m n convex hull spanned by {l m ® n; i.e. 
a// doubly stochastic matrices of the form A®B, where A G f2 m , I? G f2 n . T/ien i/ie 
extreme points of ^>m t n is the set V m <8> V n , i.e. each extreme point is of the form 
P®Q, where P G V m , Q G V n . 

Proof. Use Birkhoff 's theorem and the bilinearity of A ® B to deduce that 
Vl/ m) n is spanned by V m <8> Pn- Clearly "P m ®?„ C Pmn- Since Birkhoff 's theorem 
implies that P mn are extreme points of ft mn it follows that V m ® V n C P mn are 
convexly independent. □ 



Theorem 2.4 Let <& m ,n be the convex set of mn x mn nonnegative matrices 
characterized by 2mn+(2n— 2)m 2 + (2m— 2)n 2 linear equations of the following form. 
View C G n mnxmn as a matrix with entries C(i,k)(j,l) where i,j = l,...,m, k,l = 
1, . . . , n. Then C G H>™ nxmn belongs to <J> m , n i/ i/ie following equalities hold. 
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m,n m,n 

^ C(i,k),(j,l) = c <JM,k) = 1) » = 1, ■ ■ ■ , m, fe = 1, . . . , n, 

m m m m 

j=X j=i j=i j=i 

where i = 2, . . . , m and k, I = 1, . . . , n, 

n n n n 

1=1 1=1 1=1 1=1 

where k = 2, . . . ,n and i,j = 1, . . . ,m. 



Furthermore 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



Proof. The conditions (|2.6p state that C € Q mn . We now show the condi- 
tions ^m.n ^ <£ m .n- Let ^4 G £l m ,B G Q n and consider the Kronecker product (|2.4|) . 
Then for z, j G (m), the block of >1 ® 1? is dy-B G A n . Since A n is a cone, it 
follows that for any C E having the block form C = [C„], Cjj G M™ xn , i, j G 

(m), each Cjj G A n . Lemma 12.11 yields the conditions for each i,j G (m). Since 
A ® B = P(B <g) A)P T we also deduce the conditions ([2.7ft for each /c, Z G (n). □ 



Lemma 2.5 ^2,2 = $2,2- 

Proof. Let D = [d pq ]^ q=1 G $2,2- bmce 



^11 ■- 



it follows that 



Since 



G 



11 



it follows that 
Since 

it follows that 



-21 



C?21 ^22 



12 •= 



di3 ^14 

<^23 <^24 



G A 2 



(ill = ^22 = a, dio = d 



12 — "21 
^31 ^33 



b, = d 2 4 = c, du = & 



23 



. G12 : 



dl2 di4 
^32 ^34 



G A 2 



^31 = c, d 32 = d, C?33 = a, d 



34 



6. 



^31 CZ32 
C?4i CZ42 



22 • = 



<^33 ^34 
dis <^44 



G A 2 



C?4i = d, C?42 = C, G?43 = 6, (Z, 



11 



b. 



So a, 6, c, d > and a + 6 + c + d= 1. This set has 4 extreme points which form the 
set V 2 ®V 2 - □ 



The following result was communicated to me by J. Rosenberg. Recall that 
P G V n is called a cyclic permutation if Y2i=i P l ls a ma trix whose all entries are 
equal to 1. 
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Lemma 2.6 Let P, Q G V n be cyclic permutations. Then the block matrix D = 
^[ pi Q j ]i,j=i belongs to f n ,„. If P + Q i for i = 1, . . . , n - 1 then D * n>n . In 
particular ^f ni n S 3> n>n /or n > 4. For n = 3 eac/t .D o/ £/je above form is in ^3,3. 

Proof. Since P\Q j G Q n it follows that P l Qi G £l n for i,j = l,...,n. 
Hence the conditions (|2.8p and (|2.6p are satisfied. It is left to show the conditions 
(E2D. Denote A* = [ag]£ p=1 ,i^ = $1%=! G On- View £> = [c^y,,)]. Then 

1 ™ 

c(«,*)(j,0 = - X) a 2v > *>i»M = l,...,n. (2.10) 

P =i 

Since £)? =1 = 1 for p, Z = 1, . . . , rt and A 1 G O n we obtain \ Yuj=\ c (i,k)(j,i) = 

n ^P=l a ip = n- In a similar wa y we deduce that YTj=i c (j,k)(i,l) = So £> G $ n , n . 

Suppose that D G ^ n ,n- Observe that P n Q n = I n I n = I n . Assume D as 
a convex combination of some extreme points U ® V G V n ® P n with positive 
coefficients. Express 17® V as a block matrix [(U ® V)jj]^- =1 . Suppose furthermore 
that (C/®F) nn / nxn . Then — J n . Hence there exists j 6 (n — 1) such PQ-* — /, 
i.e P = Q n ~ 3 . If P is not a power of Q we deduce that .D G" ^n^. For n > 4 it is 
easy to construct such two permutations. For example, let P and Q are represented 
by the cycles 

l_»3_»2-»4->...->n->l, l-»2->3->4-*...-»n-»l. 

If n = 3 then one has only two cycles i? and R 2 . A straightforward calculation 
show that if P, Q G {i?, i? 2 } the D G * 3>3 . □ 

Note that the system (12. 6p has 2mn — 1 linear independent equations. Since any 
permutation matrix is an extreme point in Sl mn we deduce. 

Corollary 2.7 The convex set <& m ,n C R™ nxmn j s gwen 6y at most 2((rt - 
l)m 2 + (m — l)n 2 + ran) — 1 linear equations. It contains all the extreme points 
V m ®V n of^ m ^ n . 

It is interesting to understand the structure of the set § m ,n and to characterize 
it extreme points. It is easy to characterize the following larger set. 

Lemma 2.8 Let Q m ,n be the convex set of mn x mn nonnegative matrices 
characterized by 2mn + (2n — 2)m 2 linear equations of the following form. View 
C G W nnxmn as a matrix with entries cu,k)(j,l) where i, j = 1, . . . , m, k, I = 1, . . . , n. 
Then C G R™ nxmn belongs to 9 mi „ if the equalities Q£Q) and hold. Then 

$m,n C Q m ,n C Q mn . Furthermore, any C = [Cij]™ j=1 G @ m ,n is of the following 
form 

Cij = aijDij, Dij efi n , i, j = 1, . . . ,m, A = K,-]™ =1 G Q m . (2.11) 

In particular, the extreme points of @ m .n « r e of the the above form where A G 
V m , Dij G V n for i, j = 1, . . . , n. 
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Proof. Observe first that C in the block from C = [Cy], Cjj G M rixri where 
Cij G M" xn . Conditions (12, 8p equivalent to the assumptions that Cy G A„. Hence 
Cy = fijDij for some Djj G f2 n and > 0. If /y = we can choose any Dij G Q n . 
If > then is a unique doubly stochastic matrix. Let F = [fy] G W nxm . 
Then the conditions (|2.6p are equivalent to the condition that F G Vt m . Thus the 
conditions (12, 8j) and (|2.6h are equivalent to the statement that C = [fijDij] where 
each D^ G Q n and F = [fy] G TO . 

Since the extreme points of O n are V n we deduce that any extreme point of m>n 
is of the block form C = [fijPij] where each Py G V n . Since the extreme points 
of Q m are V m it follows that the extreme points of Q m ,n are of the form E = [Eij] 
satisfying the following conditions. There exists a permutation a : (m) — > (m) such 
that -Ejo-(j) G for i = 1, . . . , m and EJjj = mxm otherwise. □ 



3 Permutational similarity and equivalence of matrices 

For A G M. nxn denote by tr A the trace of A. Recall that (A, B), the standard inner 
product on M nxn , is given by tr AB J . 

We say that A, B G W lXn are permutationally similar, and denote it by A ~ B 
if B = PAP T . Clearly, if A ~ £> then ^4 and -B have the same characteristic 
polynomial, i.e. det(x/ n — A) = det(xl n — B). In what follows we need the following 
three lemmas. The proof of the first two straightforward and is left to the reader. 

Lemma 3.1 Let A = [a,ij],B = [bij] G M nxn . Assume A ~ B. Then the 
following conditions hold. 

P(a n , a nn ) T = (6n, . . . , b nn ) T for some P G V n , (3.1) 
R{a\2i ■ ■ • , ainj a 21, 023) • • • > 0>2n, ■ ■ ■ j«nl> • • • j«n(n-l)) = (3-2) 
(&12, • • • , bi n , &21, &23, • • • , &2n, • • • , &nl, • • • , & n.(n-l)) T / or some # G ^n 2 -n- 

Lemma 3.2 Assume that A,Be M. nxn satisfy the conditions (3.1\) and $3.2\) . 
Then tr(A + tI n )(A + tI n ) T = ti(B + tI n )(B + tI n ) T for each t G R. 

Lemma 3.3 Let A = [aij],B = [bij] G IR nxn satisfy the conditions $3. 1\) and 
113.2]) . Fix t G R such that t ^ a« — ai-k for each i,j, k G (n) such that i 7^ j. Then 
the following conditions are equivalent. 

1. A~B. 

2. B + tl n = P{A + tI n )Q T for some P,Q eV n . 

Proof. Suppose that 2 holds. Hence there exists two permutations a, n : 
(n) — > (n) such that 

bij + tSij = 0(7(i)^y) + *<*(r(i)»,Cj) for a11 J G (™)- 

Assume that a ^ rj. Then there exists i ^ j £ (n) such that cr(i) = r/(j) = fc. Hence 
bij = a/cfc + The condition (13.21) implies that fry = ai 1 j 1 for some ii 7^ j\ G (n). 
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So t = a,i 1 j 1 — akk, which contradicts the assumptions of the lemma. Hence a = rj 
which is equivalent to P = Q. Thus 

B + tI n = P(A + tl n )P T = PAP T +tl n =* B = PAP T . 
Reverse the implication in the above statement to deduce 2 from 1. □ 

We recall standard facts from linear algebra. 

Lemma 3.4 Let X = = [xi x 2 . . . x m ] G R nxm , where xi, . . . , x m G R n 

are the m columns of X. Denote by X G R mn the column vector composed of the 
columns of X, i.e. (X) T = (xj , xj, . . . , x£). Let A = [a tj ] G R mxm ,B = [b u ] G 
M. nxn . Consider the linear transformation ofW nxn to itself given by X h-> BXA t = 
[{BXAT) ki \lf =1 : 

m,n 

(BXA T ) ki = ^2 a ijbkixij, k = 1, . . . ,n, i = 1, . . . ,m. (3.3) 

Then this linear transformation is represented by the Kronecker product A®B. That 
is, 

KXA T = (A B)X for all X G R nxm . (3.4) 

Proof. Observe first that BX = [Bx± Bx2 . . . 5x m ] . This shows (|3.4p in the 
case A = I m . Consider now the case B = I n . A straightforward calculation shows 
that (A (g> L n )X = XA T . Since A® B = (A® L n )(L m ® B) we deduce the equality 
(JUj). □ 



Theorem 3.5 Let A = [a«],S = [bij] G M nxn . T/ie following conditions are 
equivalent. 

1. A~ B. 

2. The following conditions hold. 

(a) The conditions \3.1\) and \3.2\) hold. 

(b) Fix t G R such that t ^ dij — akk for each i,j,k G (n) such that i ^ j. 
Then there exists Z G ^ n .n satisfying Z(A + tl n ) = B + tL n . 

Proof. Assume 1. So B + tl n = P{A + tL n )P T for some P G V n and each 
tel. Use Lemma HQ] to deduce that (P (g> P)(A + £l n ) = B~+tl n . Hence the 
condition £6 holds. Lemma 13.11 yields the conditions (|3.ip and (|3.2|) . 

Assume £ Use Lemma O yields that tr(^ + */„)(>! + t/ n ) T = tr(B + tI n )(B + 
tL n ) T . We claim that 

^max tvP(A + tL n )Q T (B + t/ n ) T = ^max (sTIfn) T ^(-4 + t7„). (3.5) 

To find the maximum on the right-hand side it is enough to restrict the maximum 
on the right-hand side to the extreme points of ^f nn . Lemma 12.31 yields that the 
extreme points of ^ nj „ are V n ®V n - Let Y = Q ® P G V n ® V n . (|3.4p yields that 

(BTtL n ) T Y(A + tL n ) = tvP(A + tL n )Q T (B + tL n ) T . 
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Compare the above expression with the left-hand side of (|3.5p to deduce the equality 
in (|33D . 

Assume that the maximum in the left-hand side of (|3.5p is achieved for P* , Q* G 
V n . Use Cauchy-Schwarz inequality to deduce that 

trP*{A + tI n )Qj(B + tI n ) T < 
((tr P*(A + i/ n )(A + tI n ) T Pj) tr(B + i/ n )(B + tI n ) T )^ = tr(B + iJ n )(S + tI n ) T . 

Equality holds if and only if B + tl n = P*(A + tI n )Qj . The assumption 26 yields 
the opposite inequality 

tr(B + tI n ){B + tI n ) T = (B7tI n ) T Z(A + tI n ) < 
max (B + tI n ) T Y(A + tI n ) = ixP m (A + tI n )Qj (B + t/ n ) T . 

Hence B + tl n = P*(A + tI n )Qj . Lemma [3731 implies that A ~ P. □ 
The proof of the above theorem yields. 

Corollary 3.6 Assume that the conditions 2 of Theorem \3. 51 holds. Let ^ nn (^A, B) 
be the set of all Z G ^ n ,n satisfying the condition Z{A + tl n ) = B + tl n . Then all 
the extreme points of this compact convex set are of the form P®P G Vn^Vn where 
PAP T = B. 

A,B G M nxn are called permutationally equivalent, denoted as A ~ B, if B = 
PAQ T for some P G V n ,Q G Pm- The arguments of the proof of Theorem 13,51 yield . 



Theorem 3.7 Let A,B G M nxm . T/ie following conditions are equivalent. 

1. A^B. 

2. tr^ T = tvBB T and there exists Z G $ mn satisfying ZA — P. That is, 
view the entries of Z as Zu^y un where i,j G (m),k,l G (n). T/ien these 
m 2 n 2 nonnegative variables satisfy 2((n — l)m 2 + (m — l)n 2 + mn) conditions 
$2.6W2.8\) and the mn conditions: 

m,n 

Yl z (i,k)(jl) ai i = hi for k = 1,. . . ,n, i = l,...,m. (3.6) 
J\!=l 



Corollary 3.8 Assume that the conditions 2 of Theorem 3.1 holds. Let ^ m n (A, B) 



be the set of all Z G ^f m .n satisfying the condition ZA = B. Then all the extreme 
points of this compact convex set are of the form Q<S>P G V m ®V n where PAQ T = B. 



4 GIP and SGIP 

4.1 Graph isomorphisms 

Theorem 4.1 Assume that ^> n)n is characterized by fin) number of linear equal- 
ities and inequalities. Then isomorphism of two simple undirected graphs G\ = 
(V, .Ei), G*2 = (V, E2) where jfV = n is decidable in polynomial time in max(/(n), n). 
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Proof. Let A, B G {0, l} nxri be the adjacency matrices of G±,G2 respec- 
tively. Recall that A, B are symmetric and have zero diagonal. G\ and G2 are 
isomorphic if and only if A ~ B. It is left to show that the conditions 2 of Theorem 
13.51 can be verified in polynomial time in max(/(n),n). 2a means that Gi and G2 
have the same degree sequence. This requires at most An 2 computations. Assume 
that 2a holds. Note that t = 2 satisfies the first part of the condition 2b. The 
existence of Z G *& n ,n satisfying Z(A + 2I n ) = B + 2I n is equivalent to the solvabil- 
ity of f(n) + n 2 linear equations and inequalities in n 4 nonnegative variables. The 
ellipsoid method [8j [7J yields that the existence or nonexistence of such X G ^> n ,n is 
decidable in polynomial time in max(/(n),n). □ 

Theorem 4.2 Assume that ^ n>n is characterized by f{n) number of linear equal- 
ities and inequalities. Then the isomorphism of two simple directed graphs G\ = 
(V,Ei), G2 = (V,E2), (self-loops allowed), where ffV = n is decidable in polyno- 
mial time in max(/(n),n). 

Proof. Let A,B£ {0, l} nxn be the adjacency matrices of G\,G2 respectively. 
Apply part 2 of Theorem 13.51 with t = 2 to deduce the theorem. □ 

The application of part 2 of Theorem 13.51 yields. 

Theorem 4.3 Assume that ^> n . n is characterized by f{n) number of linear equal- 
ities and inequalities. Let A,B£ M. nxn . Then permutational similarity of A and B 
is decidable in polynomial time in max(/(n),n) and the entries of A and B. 

Let G = (Vi L)V2,E) be an undirected simple bipartite graph with the set of 
vertices divided to two classes V% , V2 such that E C V\ x V%. Assume that = 
n i #^2 = m and identify V\, V2 with (n), (m) respectively. Then G is represented by 
the incidence matrix A = [aij] £ {0, l} nxm where aij = 1 if and only if the vertices 
i £ {n),j G (m) are connected by an edge in E. Let H = (V\ U V2,F) be another 
bipartite graph with the incidence matrix B G {0, l} nxm . Ifm/n then G and H 
are isomorphic if and only \i Ak> B. If m = n G and H are isomorphic if and only 
if either A ~ B or A ~ B T . Theorem 13.71 yields. 

Theorem 4.4 Assume that ^ m ,n is characterized by a g(m, n) number of lin- 
ear equalities and inequalities. The isomorphism of two simple undirected bipartite 
graphs G\ = {V\ U V2, E\), G2 = (Vi U V2, E2) where jfV\ = n, V2 = m is decidable 
in polynomial time in max(g(m, n), n + m). 

Theorem 4.5 Assume that ^m^n is characterized by g(m, n) number of linear 
equalities and inequalities. Let A, B G W ixm . Then permutational equivalence of A 
and B is decidable in polynomial time in m&x(g(m, n), n + m) and the entries of A 
and B. 

We now remark that if we replace in Theorem l3.5l and Theorem 13. Ti t he sets ^n,n 
and ^ m ,n by the sets & n ,n and <& m ,n respectively, we will obtain necessary conditions 
for permutational similarity and equivalence, which can be verified in polynomial 
time. 
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4.2 Subgraph isomorphism 



Theorem 4.6 Assume that ^> ntn is characterized by f(n) number of linear equal- 
ities and inequalities. Let G3 = (W,E^),G2 = (V, -E2) be two simple undirected 
graphs, where #W = m < jfV = n. Then the problem of determining if G3 is 
isomorphic to a subgraph of G2 is decidable in polynomial time in max(/(n), n). 

Proof. Add n — m isolated vertices to G3 to obtain the graph G3 on n 
vertices. Let C, B G {0, l}™ xri be the adjacency matrices of G3, G2 respectively. We 
claim that G3 is isomorphic to a subgraph of G2 if and only if 

Z(C + 2 n2 I n ) <B + 2 n2 I n for some Z £ * n>n . (4.1) 

Assume first that G3 is isomorphic to a subgraph of G2. This is equivalent to 
the statement that PCP T < B for some P £ V n . (That is in each place where 
PCP T has entry 1, then B has entry 1 at the same place.) As PP T = I we deduce 
that flUD holds for Z = P ® P. 

Assume that ([4.ip is satisfied. Let 

Z = W ( P > Q) p ®Q> w ( p > Q) > for each PiQ^Vn and ^ w(P, Q) = 1. 

P,QeP n P,QeV n 

Hence there exists P*, Q* € V n such that w(P*,Q*) > f4w- Q4. 1 p yields that 

j^Q,(C + 2 n2 I n )Pj <B + 2 n2 I n . 

Since n = 2 n ~ for n = 1, 2 and n < 2 n ~ for 2 < n it follows that n! < 2 2 for 
n > 2. Hence (n!) 2 < 2 n for n > 1. Since all off diagonal elements of -B are at most 
1 it follows that = Q*. Hence P*CPj < (nl) 2 B. Thus if P*CPj has 1 in the 
place (i, j) then i? can not have zero in the place (i, j). That is B has 1 in the place 
Therefore G3 is isomorphic to a subgraph of G2. □ 
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